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ABSTRACT 


The  rapidly  growing  body  of  literature  on  turbine  combustor  modelling  is  reviewed. 
The  analytical  models  are  presented  and  discussed  with  particular  emphasis  placed  on  their 
ability  to  predict  gross  operating  characteristics  as  well  as  pollutant  emission  levels. 
Comparison  is  made  between  theory  and  engine  and  laboratory  experiments  showing  the 
general  inaccuracies  of  current  models.  A  suggestion  for  improvement  of  the  best  available 
model  is  made,  and  areas  for  fundamental  research  are  pointed  out. 
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SECTION  I 
INTRODUCTION 

This  survey  was  written  in  an  attempt  to  assess  the  status  of  turbine  engine  combustor 
research.  There  is  a  rapidly  growing  body  of  literature  both  on  engine  emissions  test  results 
and  on  analytical  models  of  the  combustion  process.  Prior  to  the  mid-1 960's,  the  majority 
of  research  reported  was  concerned  with  fundamental  studies  of  combustion  processes, 
and  very  little  was  published  on  the  applications  of  this  knowledge  to  practical  engines. 
As  shown  by  Gradon  (Ref.  1),  the  development  of  an  engine  relied  heavily  on  previous 
experience  and  a  great  deal  of  cut  and  weld  modifications.  At  that  time  there  was  a 
growing  emphasis  on  engine  emissions,  particularly  with  respect  to  visible  soot  and  smoke 
(Ref.  2).  The  absence  of  a  realistic  analytical  model  of  the  turbine  combustor  process 
is  adequately  emphasized  by  the  465  combustor  modifications  and  140  full-scale  engine 
tests  carried  out  on  the  Pratt  and  Whitney  JT8D  fan  engine  to  produce  a  smokeless 
combustor.  The  influence  of  the  National  Air  Pollution  Control  and  Administration 
(NAPCA)  and  the  Environmental  Protection  Agency  (EPA),  Climatic  Impact  Assessment 
Program  (C1AP)  is  being  felt  throughout  the  manufacturing  industry  and  by  users,  both 
civilian  and  military.  Research  on  automobile  emissions  has  shown  that  hydrocarbons,  nitric 
oxide,  and  sunlight  lead  to  production  of  per-oxyacetyl  nitrates  (PAN)  in  smog  (Refs. 
3  through  5).  Emission  of  carbon  dioxide,  steam,  nitric  oxides,  and  particulates  at  altitude 
has  led  to  suggestions  of  a  green-house  effect,  ozone  depletion  with  increased  cosmic  ray 
penetration  and  increased  sunlight  reflection  by  particulates  generating  a  new  ice-age. 
Accurate  assessment  of  the  real  emissions  problem  is  essential.  Required  modifications 
to  engines  should  be  based  on  predictions  of  analytical  models  coupled  with,  rather  than 
solely  by,  the  cut  and  weld  approach  if  they  are  to  be  accomplished  at  reasonable  cost. 
Currently,  analytical  models  are,  at  best,  marginally  adequate  in  predicting  emission  levels 
at  cruise  design  point  and  inadequate  in  predicting  gross  operating  characteristics  such 
as  blowout,  relight,  mixing  rates,  etc.  It  is  suggested  that  an  analytical  model  should  first 
be  capable  of  predicting  all  the  gross  characteristics  accurately  before  pollutant  emission 
predictions  are  made. 

A  typical  turbine  combustor  can  is  shown  schematically  in-  Fig.  1.  It  is  generally 
divided  into  three  zones:  primary,  secondary,  and  dilution.  The  primary  zone  consists 
of  a  fuel  injection  region  and  a  highly  stirred  reaction  region.  Combustion  occurs  at 
essentially  stoichiometric  fuel/air  ratio  (liquid  hydrocarbons  have  very  narrow  flame 
stability  limits)  with  approximately  90-percent  combustion  completion  occurring  in  it.  The 
secondary  zone  completes  the  combustion  to  about  99  percent  and  initiates  dilution, 
completed  in  the  dilution  zone,  to  provide  an  acceptable  turbine  inlet  temperature. 

The  turbine  engine,  operating  near  the  cruise  design  point,  has  a  very  efficient 
combustor,  and  overall  heat  release  rates  can  be  predicted  quite  well  (Ref.  6)  assuming 
air  is  instantaneously  mixed  with  fuel  and  burnt.  However,  performance  at  idle  (minimum 
throttle  setting,  relatively  low  turbulence  level)  and  takeoff  (maximum  throttle  setting, 
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relatively  high  turbulence  level)  are  predicted  rather  poorly  because  of  the  different 
fluidynamics  and  chemistry  interactions  in  the  combustor  at  these  throttle  settings.  A 
much  more  complicated  model  of  the  combustor  is  required,  which  considers  locally  fuel 
rich-fuel  lean  regions  caused  by  incomplete  mixing  of  fuel  and  air,  and  also  considers 
the  problems  of  blowout  and  relight.  Recent  research  has  shown  that  the  combustor  zones 
contain  several  distinct  regions  in  which  mixing  rates  are  intense  and  some  in  which  little 
mixing  occurs.  The  phenomenon  of  hysteresis  in  the  primary  zone  in  which  blowout  occurs 
at  a  higher  flow  rate  (shorter  residence  time)  than  can  be  tolerated  for  stable  ignition 
can  now  be  readily  explained.  Consider  two  high  mixing  rate  reactor  regions  coupled  by 
a  zero  mixing  rate  reactor  (plug  flow)  which  recirculates  the  products  of  one  to  the  other. 
This  additional  reactor  enables  the  high  mixing  rate  reactors  to  boot  strap  each  other 
in  such  a  way  that  the  combustion  products  of  one  helps  to  sustain  the  flame  in  the 
other.  This  model  also  indicates  that  raw  fuel  is  subjected  to  flame  temperatures  in  the 
individual  reactors  and,  in  particular,  in  the  injector  region.  Thus,  instantaneous  mixing 
and  evaporation  does  not  occur  (as  implied  in  most  analytical  models)  and  a  mechanism 
and  conditions  for  thermal  degradation  (pyrolysis)  of  fuel  to  soot  and  hydrocarbon 
fragments  obviously  exist. 


Primary 


Secondary  Dilu 
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Fig.  1  Typical  Turbojet  Combustor  Can 


Current  theoretical  models  (discussed  later)  are  inadequate  for  predicting  pollutant 
emission  levels  even  at  the  cruise  design  point.  Extensive  data  from  turbofan  engines 
(sea-level  static  tests)  were  presented  in  Ref.  7  and  correlations  of  the  data  by  Lipfert 
(Ref.  8)  have  shown  that  nitric  oxide  (NOx)  production  is  primarily  a  function  of 
combustor  inlet  temperature  and  secondarily  pressure  (see  Fig.  2).  Additonal  tests  (Ref. 
9)  with  extensive  coverage  of  the  engine  exhaust  plane  with  rakes  and  traversing  probes 
have  confirmed  these  correlations  (Fig.  3a).  A  secondary  parameter  is  specific  humidity 
as  shown  in  Figs.  3a  through  d  and  Table  I,  although  the  conclusions  of  Refs.  8  and 
9  differ  in  the  magnitude  of  its  effect.  Calculations  by  Edelman  (Ref.  10)  indicate  that 
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treating  the  primary  combustor  as  a  simple,  perfectly  stirred  reactor  (mixing  complete 
on  a  microscale  in  a  time  «  the  residence  time  -  typically  one  millisecond),  and  using 
current  NOx  chemical  reaction  schemes,  is  inadequate  to  produce  even  qualitative 
agreement  with  engine  emission  data  as  shown  in  Fig.  4.  Note  that  the  curve  was  arbitrarily 
displaced  to  coincide  with  the  match  point  and  that  the  engine  data  curve  is  concave 
upward  and  the  theory  curve  is  concave  downward.  The  qualitative  influence  of  pressure, 
temperature,  and  residence  time  on  NOx  production  is  shown  in  Fig.  5.  The  situation 
is  quite  dismal  for  prediction  of  soot  and  unburncd  hydrocarbons  for  which  no  rates 
are  known.  Lipfert  (Ref.  8)  has  provided  a  means  by  which  the  latter  can  be  estimated 
by  correlating  unburned  hydrocarbons  with  carbon  monoxide  as  shown  in  Figs.  6  and 
7.  This  particular  correlation  is,  however,  limited  to  the  class  of  engines  investigated  (for 
high  bypass,  high  pressure,  turbofans)  as  is  the  improved  version  of  Ref.  9  (see  Fig.  3d). 
A  different  correlation  can  be  expected  for  different  engines  as  evidenced  by  the  different 
correlation  for  low  pressure,  turbojet  engines  (Refs.  11  and  12). 
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Fig.  2  Nitric  Oxide  Emission  Data  Correlation  Corrected  to  0.01  Humidity 


3 


AEDC-TR-73-163 


a.  JT3D  Environmental  Test— Oxides  of  Nitrogen  versus  Burner 
Inlet  Temperature  {Corrected  Data,  High  Power) 


0  20  40  60  80  100 

Engine  Inlet  Temperature,  Tjj.  °F 

b.  Correction  Factor  Curves  for  CO 
Fig.  3  JT3D  Emission  Data  Corrected  for  Humidity 
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TABLE  I 

EMPIRICAL  FORMULAS  FOR  EMISSION  LEVELS 


Oxides  ol  Nitrogen 

EIn0  ■  0.481  +  0.8?8x  10"3Fn  -  0.3655  x  10'7  FN2 

+  0  1513  x  1CT 11  Fn3  +  0. 03371(2  -  101. 18  SP  HMO 

Carbon  Monoxide 

(1/EIC0!1/5  ■  0.7148  ♦  0  355  x  10‘4FN  +  0.177  x  10‘7T|? 

+  3.797  SP  HMD 

Total  Hydrocarbons 

«1/E|thc,1,I°  ■  °-5481  +  0.360  x  10‘4  Fn  +  9  699  SP  HMO 

Where  El  ■  Emissions  Index.  IW1000  lb  Fuel 
Ffl  •  Engine  Thrust,  lb 

T(j  ■  Ambient  llnleti  Temperature,  °F 

SP  HMD  ■  Specific  Humidity,  IbH^OflbAir 
The  models  were  also  used  to  develop  correction  factor 
curves  (or  humidity  and  temperature  changes  where 
appropriate.  These  are  shown  in  Figs.  1,  and  3a 
through  d. 

(High  Pressure  Enginesl 


Fig.  4  Composite  Nitric  Oxide  Emission  Index 
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Fig.  5  Concentration  of  Nitric  Oxide  as  a  Function  of  Residence  Time  for 
Various  Temperatures  and  Pressures— Adiabatic  Reactor 


Carbon  Monoxide  Concentration,  Cco.  ppm  (Ref.  B) 

Fig.  6  Hydrocarbon/Carbon  Monoxide  Emission  Correlation 
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Burner  Inefficiency  (1  -  ijg) 


Fig.  7  Prediction  of  Carbon  Monoxide  and  Hydrocarbon 
Emission  Levels  from  Burner  Efficiency 

The  importance  of  a  good  theoretical  model  which  can  account  for  the  nonlinear 
response  of  combustor  processes  to  changes  in  design  and  operating  conditions  is 
emphasized  in  the  following  example:  "Smokeless"  combustors  were  designed  primarily 
by  promoting  better  fuel/air  mixing  (injecting  primary  zone  air  into  the  fuel  sprays  and 
improving  the  distribution  of  liner  cooling  air  to  minimize  hot  spots,  etc.).  This  reduced 
fuel  rich  regions;  however,  the  higher  average  flame  temperature  expected  was  assumed 
to  increase  NOx  production.  Lipfert  (Ref.  8)  showed  that  when  sufficient  data  were 
statistically  averaged,  the  expected  increase  was  not  apparent.  It  will  be  shown  later  that 
there  is  an  almost  complete  inability  to  predict  NOx  levels  under  fuel  rich  conditions; 
the  observed  effects  of  humidity  are  also  not  predictable.  Thus,  it  can  be  concluded  that 
current  models  are  a  poor  representation  of  the  details  of  most  features  of  turbine 
combustor  operation. 

The  following  sections  outline  some  of  the  research  that  has  been  conducted  and 
attempts  to  define  what  combination  of  different  approaches  might  lead  to  more 
quantitative  accuracy  in  predicting  pollutant  emission  levels  and  turbine  combustor 
performance  over  the  range  of  throttle  settings.  It  is  suggested  that  a  model  which  predicts 
observed  trends  over  the  whole  performance  range  is  the  most  urgent  task  to  be  solved 
before  realistic  design  changes  may  be  recommended  (to  modify  pollutant  emission  levels) 
to  the  machine  shop. 
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The  primary  problem  in  modelling  any  practical  combustion  system  is  to  account 
for  turbulent  mixing  and  chemical  reaction  simultaneously.  Although  mixing  models  have 
received  intense  interest  over  the  years  (Harsha,  Ref.  13),  it  is  only  recently  that  an 
adequate  representation  of  several  classes  of  mixing  flows  (without  combustion)  could 
be  obtained  from  one  theory  (Refs.  14  through  18).  However,  notable  as  these  successes 
are,  the  theory  is  only  applicable  to  simple  flows  involving  two  dimensions  (planar  or 
axial  symmetry).  Recently,  progress  has  been  made  in  techniques  for  including  finite  rate 
chemistry  in  a  turbulent  combustion  model  (again,  a  simple  flow  configuration)  in  which 
the  time-average  density  in  a  nonhomogeneous  reacting  flow  is  obtained  through  an 
averaging  model  rather  than  as  a  function  of  the  time-average  pressure,  temperature,  and 
species  concentration  (Rhodes  and  Harsha,  Ref.  19). 

Combustors  involve  three-dimensional  recirculating  flows.  Currently,  no  analytical 
model  is  capable  of  accurately  representing  such  three-dimensional  flows  with  recirculation. 

A  technique  for  two-dimensional  flows  developed  by  Gosman,  Ref.  20,  is  being  actively 
developed  by  many  research  groups  but  is  several  years  away  from  being  a  design  tool 
for  other  than  simple  flows. 

Currently,  combustors  are  modelled  by  three  groups  of  theories.  The  first  group  is 
based  on  calculating  an  augmented  laminar  flame  speed  (turbulent  flame  speed)  as  a 
function  of  the  turbulence  parameters  of  the  flow.  The  chemical  reaction  rate  is  introduced 
by  the  use  of  this  augmented  laminar  flame  speed  in  laminar  flame  theory.  The  second 
group  develops  a  microvolume  burning  concept  by  identifying  eddies  of  size  and  lifetime 
dependent  on  the  turbulence  parameters.  These  eddies  are  assumed  to  bum  either  from 
the  surface  at  the  laminar  flame  speed  or  homogeneously  throughout  the  volume.  The 
third  group  of  models  assumes  that  the  turbulence  level  is  so  high  that  the  mixing  time' 
can  be  either  neglected  or  is  very  small  compared  to  the  chemical  reaction  time  (perfectly 
stirred  reactor  models).  The  combustion  system  is  usually  modelled  as  a  combination  of 
stirred  reactor  volumes.  In  a  recent  analysis  (Ref.  21)  the  mixing  rate  is  considered  to  ' 
be  fast  but  finite  compared  to  the  chemical  reaction  rate  (well-stirred  reactors)  and,  in 
addition,  includes  plug  flow  in  the  primary  zone  to  preserve  the  known  hysteresis 
characteristics  associated  with  blowout  and  relight. 

2.1  TURBULENT  FLAME  SPEED  MODELS 

Turbulent  flame  speed  models,  developed  primarily  in  Russia  (Refs.  22  and  23) 
account  for  the  effect  of  turbulence  by  calculating  a  turbulent  flame  speed  based  on  analogy 
with  the  laminar  flame  speed  (Refs.  24  and  25).  This  turbulent  flame  speed  is  then 
employed  to  calculate  flame  propagation  from  various  ignition  sources  for  determining 
required  overall  combustor  dimensions  and  performance.  The  analysis  has  been  developed 
for  ramjet-type  combustors  of  constant  cross-sectional  area  with  premixed  combustibles 
and,  hence,  is  not  strictly  applicable  to  turbine  combustors. 
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In  the  highly  turbulent  combustor  situation,  burning  is  postulated  to  occur  in  a 
discrete  zone.  The  forward  boundary  of  this  zone  corresponds  to  ignition  and  is  called 
the  "flame  front."  The  aft  boundary  of  the  zone  is  marked  by  the  cessation  of  heat 
release;  therefore,  the  zone  thickness  is  determined  by  the  chemical  reaction  rate  (implying 
streamline  flow  through  the  zone)  and  gas  transport  properties.  Combustor  design  is  based 
on  determining  the  size  of  the  zone  and  its  extent  relative  to  the  flameholder.  The  major 
deficiencies  of  this  model,  according  to  Hammond  (Ref.  26).  are  the  assumptions  of 
smooth,  premixed  flames  instead  of  the  more  realistic,  highly  irregular  turbulent  flames 
in  which  mixing  plays  a  significant  role  as  in  practical  turbine  combustors.  In  addition, 
the  turbulent  flame  speed  has  not  been  established  as  a  proper  characterization  of  the 
combustion  process,  and  the  exact  dependence  of  the  turbulent  flame  speed  on  reactant 
properties  and  combustor  flow  conditions  has  not  been  established. 

2.2  MICROBURNING  MODELS 

The  microvolume  burning  concept  was  introduced  by  Shchetinkov  (Ref.  27)  and 
represents  an  attempt  to  model  both  turbulent  mixing  and  burning  simultaneously. 

Shchetinkov  postulated  mixing  to  be  characterized  by  the  formation  and  dissipation 
of  turbulent  eddies.  These  eddies  are  reasonably  discrete  masses  of  fluid  which  move  as 
a  particle  relative  to  other  eddies.  Eddies  containing  burned  or  partially  burned  gases  are 
carried  out  of  the  combustion  zone  by  turbulent  fluctuations,  after  which  they  mix  with 
eddies  of  fresh  gases.  If  the  resulting  eddies  possess  a  sufficiently  high  temperature,  ignition 
will  occur  and  the  eddy  will  burn.  Shchetinkov  modelled  the  mixing  behavior  by  means 
of  a  grid  of  microvolume  cubes  of  side  22  where  £  is  the  mixing  length.  The  lifetime 
of  the  eddies  was  defined  as  t  =  fi/u  where  u  is  the  turbulent  intensity.  The  time  available 
for  combustion  is 


tc  t  *  tjn  -  tj 

where  tm  is  a  measure  of  the  mixing  rate  within  the  eddy  and  t;  is  an  ignition  delay 
time.  Within  the  eddy,  diffusion  processes  operate  to  equalize  chemical  composition  and 
temperature;  however,  homogeneity  is  not  always  assured  and  both  surface  and  volume 
combustion  can  occur.  Shchetinkov  considered  only  volume  combustion  (homogeneous) 
while  Berl,  el  al.  (Ref.  28)  treated  both. 

This  type  of  theory  approximates  the  physical  system  rather  closely,  since  in  turbulent 
flow,  eddies  are  formed  which  travel  a  certain  mixing  length  and  decay.  However,  uniformly 
sized  eddies  rather  than  a  statistical  size  distribution  is  a  poor  assumption  and  theories 
for  calculating  mixing  lengths  are  poorly  developed.  A  more  promising  approach  to  the 
details  of  combustion  in  turbulence  is  that  of  Rhodes  and  Harsha  (Ref.  1 9).  They  consider 
eddy  classes  to  be  formed  with  each  eddy  class  corresponding  to  a  prescribed  concentration 
and  a  probability  of  existence  at  each  point  in  the  flow.  The  probability  distribution 
functional  form  is  assumed  to  be  known;  thus  the  mean  and  standard  deviation  define 
the  value  of  the  function.  Assuming  a  Crocco  relation  between  concentration  and  velocity, 
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similar  to  that  between  velocity  and  enthalpy,  the  standard  deviation  is  defined  in  terms 
of  the  axial  velocity  fluctuation  which  is  related  to  the  turbulent  kinetic  energy.  Chemical 
reaction  is  incorporated  by  considering  concentration  classes  as  transient,  perfectly  stirred 
reactors  (TPSR).  A  TPSR  is  one  in  which  the  feed  rate  and  feed  concentration  may  vary 
with  time,  and  the  time  scale  considered  is  short  compared  to  the  time  required  for  the 
reactor  to  reach  a  steady-state  condition.  However,  much  more  work  is  required  to  make 
this  a  viable  theory  for  flows  other  than  the  simple  axisymmctric  flow  considered  in  Ref. 
19. 

2.3  STIRRED  REACTOR  MODELS 

A  perfectly  stirred  reactor  (PSR)  is  a  region  containing  a  homogeneous  reacting 
mixture;  therefore,  if  samples  were  taken  from  arbitrary  locations  in  the  zone,  they  would 
be  identical  in  temperature,  pressure,  and  chemical  composition.  The  feed  material  is 
instantaneously  mixed  throughout  the  volume,  thus  maintaining  the  uniform  composition 
of  the  zone,  and  an  equal  mass  flow  of  exhaust  products  leaves  the  zoneliavihg  the  same 
composition  as  the  reactor  contents,  The  performance  of  the  combustor  is  determined 
solely  by  chemical  kinetics  and,  in  general,  the  reactions  do  not  progress  to  completion 
in  the  reactor. 

The  PSR  concept  was  first  applied  to  combustor  modelling  by  Longwell  and  Weiss 
(Ref.  29).  Initial  attempts  were  based  on  an  assumed  single  PSR  for  the  combustor  (Refs. 
30  through  35).  However,  in  a  practical  combustor  (discussed  later),  the  physical 
phenomena  of  blowout,  a  residual  flame,  hysteresis  between  blowout  and  relight,  flame 
instability,  incomplete  mixing,  etc.,  require  more  sophisticated  models.  Qualitative 
predictions  of  some  of  these  phenomena  were  obtained  by  several  authors  (Refs.  34  through 
38).  primarily  with  the  inclusion  of  recirculation  of  reaction  products  to  one  or  more 
stirred  reactors.  The  inclusion  of  detailed  finite  rate  chemistry  to  a  single  PSR  by  Spalding 
Ref.  39,  Jones  (Ref.  40);  and  Jenkins  (Ref. 

41)  provided  the  basis  for  the  prediction  of 
pollutants  from  turbine  engines.  Hammond 
and  Mellor  (Ref.  26)  presented  results  of 
calculations  assuming  the  turbine  combustor 
could  be  modelled  by  a  primary  combustor 
(Fig.  8)  consisting  of  two  PSR's  in  parallel, 
i.e.,  a  second  reactor  recycles  part  of  the 
output  of  the  first  reactor;  a  secondary 
combustor,  consisting  of  a  series  of  plug-flow 
reactors  (PFR)  each  fed  with  secondary  air 
instantaneously  mixed  with  the  output  of  the 
previous  reactor,  followed  by  a  final  PFR  fed 
with  dilution  air.  As  usual,  the  reactor 
volumes  and  recirculation  rates  were  chosen 
arbitrarily.  The  injected  fuel  is  all  delivered 
to  the  first  reactor  with  burnout  achieved  in 
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Fig.  8  Hammond's  Recirculation 
Reactor  Model 
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the  PFR's  in  the  secondary  combustor.  The  model  was  modified  by  assuming  part  of 
the  fuel  was  delivered  to  the  second  reactor  (Ref.  42).  Neither  model  is  capable  of 
predicting  a  characteristic  delay  time  obtained  when  a  recirculating  PFR  is  included  in 
the  model,  observed  by  Beer  and  Lee  (Ref.  43)  and  shown  to  be  essential  for  predicting 
flame  stabilization  phenomena  (Refs.  26,  37,  21,  and  44).  This  delay  time  is  observed 
in  transient  flow  tests  in  which  a  trace  gas  in  the  feed  is  suddenly  cut  off.  Thus  in  a 
plug-flow  recirculating  model,  tracer  material  is  continuously  recycled  back  through  the 
stirred  reactors.  A  trace  gas  monitor  shows  a  characteristic  delay  time  before  tracer  gas 
concentration  follows  the  1/e  decay  of  stirred  reactors. 

The  work  reported  by  Swithenbank  and  co-workers  (Refs.  21  and  44)  represents 
the  most  complete  analysis  to  date,  with  the  exception  of  inability  to  calculate  pollutant 
concentration.  Their  model  is  a  development  of  the  basic  theory  of  Vulis  (Ref.  35)  and 
Swithenbank  (Ref.  45)  in  which  the  primary  combustor  is  modelled  by  two  well-stirred 
reactors  (WSR)  in  series.  The  well-stirred  reactor  is  one  in  which  the  fuel  and  oxidant 
are  mixed  rapidly  but  not  instantaneously.  The  two  WSR's  are  coupled  by  a  PFR  which 
recirculates  part  of  the  second  stirred  reactor  to  the  feed  stream  of  the  first  (Fig.  9). 
The  stirred  reactors  are  modelled  as  well-stirred  reactors  by  incorporating  a  mixing  delay 
time.  The  mixing  time  is  obtained  from  considerations  of  the  turbulence  level  generated 
by  total  pressure  losses  in  the  air  and  fuel  injection  processes.  In  addition,  entrainment 
theory  is  used  to  determine  estimates  of  the  individual  reactor  volumes  and  mass  flows. 
More  fundamental  well-stirred  reactor  theories  have  previously  been  developed  (Refs.  46 
through  52),  but  none  of  these  have  been  incorporated  into  more  than  a  single  reactor. 
They  may  be  a  logical  extension  of  the  Swithenbank  model  with  the  latter  used  for 
predicting  the  turbulence  level,  reactor  volumes,  mass  flows,  and  hence  residence  times 
in  a  complex  combustor  such  as  in  the  turbojet.  This  model  is  readily  extended  to  include 

* - * - - 
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Fig.  9  Swithenbank's  Recirculation  Reactor  Model 
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detailed  pollutant  concentrations  since  the  necessary  data,  reactor  volumes,  and  volume 
flow  rates  (hence  reactor  residence  times)  are  determined.  The  utility  of  this  model  is 
considered  sufficiently  promising  that  details  of  the  calculation  of  reactor  volumes,  flow 
rates,  and  some  experimental  evidence  to  substantiate  the  physical  model  are  included. 
The  ability  to  calculate  individual  reactor  residence  time  is  unique  to  this  model,  and 
the  importance  of  this  variable  is  shown  in  Fig.  5. 


2.4  SWI  THEN  BANK  MODEL 


The  Swithenbank  modeKis  the  result  of  several  years  of  experience  in  modelling 
practical  combustors  and  is  oriented  at  a  specific  model  of  a  turbojet-type  combustor. 

The  experimental  combustor,  used  by  Swithenbank,  is  shown  in  Fig.  10,  and  the 
reactor  volumes  are  shown  schematically  in  Fig.  11.  Rather  than  arbitrarily  choosing- the 
reactor  volumes  as  did'  Hammond  (Ref.  42),  the  method  of  Swithenbank  is  outlined  to 
illustrate  how  they  may  be  calculated. 


1.00 


Fuel 


Mass  Flow  Rate. 


Fig.  11  Lycoming  2.75-in.  Can:  Proposed  Interconnection  of  Reactors 
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2.4.1  Primary  Zone 

The  primary  combustor  is  modelled  as  a  pair  of  stirred  reactors  coupled  by  a  plug-flow 
reactor.  The  volumes  and  mass  flows  of  the  reactors  are  calculated  using  jet  entrainment 
theory.  The  total  primary  zone  volume  is  obtained  by  experiment  or  by  inspection  of 
the  actual  combustor  being  investigated.  Swithenbank  used  a  salt-tipped  probe  to  determine 
the  position  on  the  combustor  axis  at  which  the  whole  flow  field  (primary  zone)  became 
sodium  colored.  If  this  procedure  is  not  possible,  the  midpoint  between  the  primary  and 
secondary  air  inlet  holes  may  be  chosen.  Although  the  swirler  contributes  to  the  primary 
flow,  the  primary  air  holes  are  largely  responsible  for  the  formation  and  intensity  of  the 
recirculation.  The  amount  of  recirculation  can  be  predicted  quite  simply  as  follows: 

The  amount  of  flow  entrained  into  a  jet  is  given  by  (Ref.  53) 


dj  =  d0  Cd 


where  liij  is  the  mass  flow  of  air  in  the  jet  entering  via  an  orifice  of  diameter  d„  and 
discharge  coefficient  Cp  ,  Pj  is  the  air  density,  and  x  is  the  distance  penetrated  by  the 
jet  into  the  combustor.  The  distance  available  to  entrain  flow  is  the  radius  of  the  can. 

Swithenbank  assumed  pjp  =  4;  however,  p  could  readily  be  calculated  assuming 
equilibrium  combustion  of  90  percent  of  the  fuel  in  the  primary  zone.  From  Fig.  10, 
x  =  1.375  in.,  d0  =  0.26  in.,  and  Cp  is  given  as  0.85;  thus 


■£r  =  °-92 

mj 


The  primary  region,  therefore,  behaves  as  a  stirred  reactor  with  92  percent  of  the 
primary  air  mass  flow  entrained.  Only  a  fraction  (7.8/25.5  =  0.3)  of  this  required 
entrainment  mass  flow  is  available  from  the  swirler;  thus  62  percent  of  the  primary  flow 
must  be  recycled  into  the  swirler  region.  Measurements  and  water  models  of  this  type 
of  flow  have  shown  it  to  be  relatively  smooth  flow:  hence,  it  is  treated  as  a  plug-flow 
reactor.  Addition  of  recycle  flow  to  swirler  flow  takes  place  in  the  ratio  0.62/0.3  =  2.06. 
The  strong  shearing  region  in  the  swirler  flow  is  treated  as  a  stirred  reactor. 


The  volumes  of  the  individual  reactors  is  estimated,  assuming  the  gas  density  is  the 
same  in  all  the  reactors:  primary,  recycle,  and  swirler.  The  volumes  of  the  swirler  and 
recirculation  reactors  are  ratioed  according  to  the  mass  flows  in  the  reactors  and  that 
of  the  primary  assuming  it  is  symmetrical  about  the  primary  orifices. 
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2.4.2  Secondary  and  Dilution  Zones 

Ihe  secondary  stirred  reactor  is  also  assumed  to  be  symmetrical  about  the  secondary 
orifice  with  its  upstream  boundary  abutting  the  primary  zone.  The  dilution  zone  stirred 
reactor  is  assumed  to  have  the  same  volume  for  equal  air  mass  flow.  The  secondary  plug 
flow  reactor  lies  between  the  secondary  and  dilution  stirred  reactors.  There  is  no  upstream 
recirculation  in  the  secondary  zone  because  there  is  more  than  sufficient  mass  flow  to 
satisfy  entrainment  requirements: 


Swithenbank  assumes  p-Jp  =  5  in  the  secondary  zone,  giving  m/mj  =  0.76  i.e.,  the  entrained 
flow  is  76  percent  of  the  secondary  airflow.  The  mass  flow  upstream  is  (25.5  +  7.8)/29.9 
=  1.08  times  the  secondary  airflow  and,  hence,  is  more  than  sufficient  to  satisfy  the 
entrainment  appetite.  The  jet  turbulence  is  assumed  to  dissipate  rapidly  (within  a  duct 
diameter)  because  of  impingement  of  radially  opposed  jets;  thus,  little  additional  mixing 
should  occur  downstream  of  the  secondary  stirred  reactor  and  the  flow  is  assumed  to 
be  plug  flow.  The  incomplete  mass  flow  entrainment  should  result  in  a  dip  in  the 
temperature  profile  at  the  axis.  The  dilution  zone  air  entrains  about  77  percent  (pj/p 
=  4)  with  1.72  times  the  dilution  air  mass  flow  upstream.  A  significant  dip  in  the 
temperature  profile  should  be  observed.  Temperature  measurements  by  Swithenbank 
confirmed  this  cold  core  as  shown  in  Fig.  12. 
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2.4.3  Fully  Mixed  Calculations 

The  mass  flows,  volumes,  and  residence  times  (volume/volume  flow  rate)  having  been 
determined,  complete  chemical  kinetic  calculations  can  be  made.  Additional  input  data 
required  are  inlet  pressures  and  temperatures.  The  calculation  can  be  initiated  assuming 
equilibrium  stoichiometric  products  from  the  second  stirred  reactor  and  plug-flow  reactor 
in  the  primary  zone  are  mixed  with  fuel  and  air  from  the  swirler  as  feed  to  the  first. 
The  output  calculated  for  the  first  reactor  is  fed  to  the  second  reactor  together  with 
the  primary  airflow.  Part  of  the  flow  output  is  then  recycled  through  a  plug-flow 
calculation.  The  whole  process  is  repeated  until  no  further  composition  changes  occur 
in  each  reactor.  Calculations  are  then  resumed  for  the  secondary  and  dilution  zone  reactors. 
A  comprehensive  set  of  reactions  and  reaction  rates  is  given  in  Table  II. 

TABLE  II 

CHEMICAL  REACTIONS  IN  THE  C/H/N/O  SYSTEM 


Reaction 

Reaction  Rate 

A 

n 

E 

Ref. 

1 

H  +  C02  ■  CO  +  OH 

5.6  x  1013  -. 

0 

23.500 

106 

2 

0  ♦  N2  -  NO  ♦  N 

1.44  x  1014 

0 

75.580 

107 

3 

0  ♦  NO  -  N  +  Oj 

4. 1  x  10* 

1 

38.340 

107 

4 

N  ♦  OH  ■  NO  ♦  H 

4.21  xlO13 

0 

0 

108 

5 

OH  +  H2  -  H^O  +  H 

2.2  x  1013 

0 

5,150 

106 

6 

OH  ♦  OH  -  0  +  HpO 

5. 75  x  lO17 

0 

780 

106 

7 

O  +  ty-  H  +  0H 

1.74  xlO13 

0 

9.450 

106 

8 

H  +  Oj  ■  0  ♦  OH 

1.56  x  1014 

0 

16.633 

109 

9 

O  +  H  +  M-OH  +  M 

3.6  x  1017 

-1 

0 

no 

10 

o  +  o  +  m-o2  +  m 

3.6  x  1017 

-1 

0 

109 

11 

h  +  h  +  m-h2*m 

1.81  x  101* 

-1 

0 

109 

12 

H  ♦  OH  +  M  •  HjO  +  M 

7.3  x  101* 

-1 

109 

13 

0  +  C02  -  CO  +  o2 

1.9  x  10*3 

0 

54. 150 

106 

14 

CO  +  0  +  M  ■  C02  ♦  M 

6.0x  1017 

-1 

2.484 

110 

15 

N2  +  Oj  -  N  +  N02 

2.7  xlO14 

■1 

120.428 

56 

16 

N2  +  02  ■  NO  +  NO 

4.2  x  1074 

0 

119, 100 

106 

17 

NO  +  NO  •  N  +  N02 

3. Ox  1011 

0 

0 

106 

18 

NO  ♦  M  •  0  ♦  N  +  M 

2.27 x  1017 

-0.5 

148.846 

56 

19 

0  +  NO  +  M  •  N02  ♦  M 

1.05  x  1015 

0 

-1,870 

105 

20 

N  ♦  02  +  M  ■  NOj  +  M 

7.0  x  1011 

-1 

56 

21 

0  +  N02  ■  NO  +  02 

1.0  x  1013 

0 

600 

106 

22 

H  ♦  N02  -  NO  ♦  OH 

7. 25  x  1014 

0 

1.930 

• 

23 

N  ♦  C02  -  CO  +  NO 

2.0x  1011 

-0.5 

7.950 

• 

24 

CO  +  N02  ■  NO  *  C02 

2. Ox  I011 

-0.5 

4.968 

• 

25 

H  +  N^  -  OH  ♦  N2 

3.01  x  1013 

0 

101800 

108 

26 

0  ♦  •  02  ♦  N2 

3.61  *  1013 

0 

24,800 

107 

27 

N2  ♦  N02  ■  NO  +  NjO 

0 

83,000 

107 

Reaction  Rate  Constant  k  -  ATn  exp  f-E/RII 

Reverse  Reaction  Rate  Constant  Obtained  from  k  and  the  Equilibrium  Constant 
•Recommended  by  R  B.  Edelman 
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TABLE  II  (Concluded) 


Reaction 

Reaction  Rite 

A 

n 

£ 

Ref. 

28 

NO  +  HOj  ■  OH  ♦  NO2 

6.0*  10*1 

0 

0 

106 

29 

0  +  N?  *  M  •  NjO  ♦  M 

6.3*  1014 

0 

56.800 

107 

30 

H  ♦  Oj  +  M  ■  H02  ♦  M 

1.5*  1015 

0 

•1,000 

106 

31 

H  ♦  Ht>2  •  OH  ♦  OH 

2.5  *  I014 

0 

1.900 

106 

32 

OH  +  HOj  ■  HjO  ♦  Oj 

1.2*  1013 

0 

1.000 

106 

33 

0  ♦  H0Z  ■  OH  ♦  Oj 

5.0*  1013 

0 

1.000 

106 

34 

H  ♦  H02  ■  Hj  +  Oz 

2.5  *  1033 

0 

700 

106 

35 

H  ♦  HOj  •  0  ♦  HjO 

1.0*  1013 

0 

1.000 

106 

36 

H7  ♦  HOj  ■  H  + 

1.9*  10° 

0 

24.000 

106 

37 

Hj02  +  M  -20H  +  M 

7.1  *  1014 

0 

-5.100 

106 

38 

Cfy  ♦  +  H  •*  M 

1.5*  ID19 

0 

99,960 

72 

39 

CH4  +  H  •  CHj  +  Hz 

5.1*  1013 

0 

12.900 

109 

40 

CH4  ♦  OH  ■  CHj  ♦  H2O 

2.85*  ID13 

0 

4.968 

72 

41 

CH4  +  0  -  CHj  +  OH 

1.7  x  ID13 

.0 

8,700 

72 

42 

CH3  +  0;  •  HCO  ♦  HjO 

1.0*  1011 

0 

0 

72 

43 

HCO  +  OH  ■  CO  ♦  H^O 

3.0*10° 

0 

0 

111 

44 

H  ♦  CO  +  M  ■  HCO  +  M 

1.0*  1017 

-1 

0 

111 

45 

HOj  ♦  HOj  -  H^  ♦  Oj 

6.5*  1013 

0 

0 

106 

46 

HjO;  +  H  ■  HpO  +  OH 

3.18*  1014 

9.000 

106 

47 

CH2*  1/202 -CO*  H2+ 

5.52*  10*p‘OE5 

1 

24.642 

79 

48 

H  ♦  NO  +  M  ■  HNO  ♦  M 

4.8*  10° 

0 

112 

49 

HNO  ♦  H  ♦  M  -  HjNO  ♦  M 

50 

HNO  ♦  HNO  ♦  M  -  H2O  +  N20  +  M 

51 

R  ♦  NO  +  M  -  RNO  ♦  M 

52 

R  +  RNO^M-RpNO  +  M 

53 

R  +  R2NO  ♦  M  -  R^NOR  +  M 

1.4*  1011 

54 

H  *  HNO  -  Hz  +  NO 

0.5 

113 

55 

H  ♦  HCO  ■  H2  +  CO 

2.0*  1013 

0 

0 

111 

+Any  CnHm  hydrocarbon  is  given  this  rate,  however,  with  n  CO  molecules  and  m/2  H? 
molecules  produced  at  this  rate. 

2.5  INCOMPLETE  MIXING  CALCULATIONS 

At  this  point  it  is  not  clear  how  to  account  for  incomplete  mixing.  Both  Heywood 
(Ref.  54)  and  Bowman  (Ref.  52)  have  carried  out  calculations,  but  Bowman  has  indicated 
deficiencies  in  his  method.  Heywood  has  not  carried  out  calculations  with  more  than  one 
reactor  in  the  primary  zone;  thus  it  is  not  clear  how  more  reactors  would  influence  the 
results.  The  method  of  Rhodes  (Ref.  19)  requires  some  estimate  of  turbulence  intensity 
and  Swithenbank's  (Ref.  43)  pressure  loss  calculation  would  supply  an  estimate.  The  rate 
of  decay  of  turbulence  could  be  obtained  by  the  method  of  Corrsin  (Ref.  46)  as  a  first 
approximation. 

Results  of  calculations  with  the  fully  mixed  model  should  be  generated,  followed 
by  an  attempt  to  include  the  above  incomplete  mixing  approaches  to  individual  reactor 
volumes.  Comparison  of  these  results  and  with  experimental  data  should  determine  what, 
if  any,  improvement  in  prediction  capability  has  been  achieved. 


17 


AEDC -TR-73-163 


SECTION  III 

POLLUTANT  CHEMISTRY  MODELS 


Several  practical  mathematical  techniques  now  exist  for  solving  the  nonlinear  algebraic 
equations  which  describe  stirred  reactor  combustion  for  complete  chemical  kinetic  schemes 
(Refs.  52,  55,  and  56).  A  large  number  of  solutions  of  problems  involving  pollutant 
production  have  been  published  (Refs.  26,  52,  and  56  through  71).  The  analyses  are 
characterized  by  their  diverse  opinions  on  the  chemical  reactions  considered  sufficient  to 
describe  the  combustion  process.  Some  reactions  have  been  included  in  all  analyses, 
although  some  differences  in  reaction  rate  constants  are  evident.  These  reactions  include: 


CO  Oxidation: 


1 

CO  • 

H  OH 

=  H  +  C02 

NO 

Formation: 

_(Zel'dovich  mechanism) 

2 

0  + 

N2  = 

NO  +  N 

3 

N  + 

02  = 

NO  +  0 

4 

N  + 

OH  = 

*  NO  +  H 

H? 

Oxidation: 

5 

OH 

+■  h2 

=  H20  +  H 

6 

OH 

+  OH 

-  o  +  h2o 

7 

0  + 

h2  = 

H  +  OH 

8 

H  + 

02  = 

0  +  OH 

9 

0  + 

H  + 

M  *  OH  +  M 

10 

0  + 

0  + 

M  =  02  +  M 

11 

H  + 

H  + 

M  =  H2  +  M 

12  H  +  OH  +  M  =  H20  +  M 


Additional  reactions  which  have  been  included  in  various  analyses  are: 


13 

CO  t  02  =  C02  +  0 

14 

CO  +  0  +  M  —  C02  + 

M 

15 

N2  +  02  =  N  +  N02 

16 

N2  +  02  =  NO  +  NO 

17 

N  +  N02  =  NO  +  NO 

18 

0  +  N  +  M  =  NO  +  M 

19 

0  +  NO  +  M  =  N02  + 

M 

18 


AEDC-TR-73-1  63 


20 

N  +  02  + 

M  =  N02  +  M 

21 

NO  +  02  : 

=  0  +  no2 

22 

H  +  N02 

=  NO  +  OH 

23 

N  +  C02  = 

=  CO  +  NO 

24 

CO  +  no2 

=  NO  +  C02 

Additional  reactions  which  have  been  considered  but  not  included  in  published 
calculations  are: 


25 

H  +  N20  = 

=  OH 

+  n2 

26 

O  +  N20  = 

=  o2 

+  n2 

27 

N2  +  no2 

=  NO  +  N20 

28 

NO  +  H02 

=  OH  +  N02 

29 

O  +  n2  + 

M  = 

N20  +  M 

30 

H  +  02  + 

M  = 

H02  +  M 

31 

H  +  H02  = 

=  OH 

+  OH 

32 

OH  +  H02 

-  H20  +  02 

33 

O  +  ho2  = 

=  OH 

04 

O 

+ 

34 

H  +  H02  = 

=  H2 

+  o2 

35 

H  +  H02  = 

=  O  H 

H  H20 

36 

H2  +  H02 

=  H 

+  H202 

37 

H202  +  M 

=  OH  +  OH  +  M 

3.1  HYDROCARBON  COMBUSTION 

A  number  of  hydrocarbons  have  been  investigated  from  methane  to  JP4  and  JP5 
fuels.  The  mechanism  for  methane  combustion  has  received  the  most  attention  because 
of  its  simplicity  compared  to  the  higher  molecular  weight  hydrocarbons. 


38 

ch4 

+  M  =  CH3  +  H  +  M 

39 

ch4 

+  H  =  CH3  +  H2 

40 

ch4 

+  OH  =  CH3  +  H20 

41 

ch4 

+  0  =  CH3  +  OH 

42 

ch3 

+  02  =  HCO  +  H20 

43 

HCO  +  OH  =  CO  +  H20 

44 

HCO 

+  M  =  H  +  CO  +  M 

This  mechanism  was  suggested  by  Marteney  (Ref.  72)  together  with  reactions  1,  2, 
3,  18,  5,  6,  7,  and  8.  The  mechanism  was  expanded  by  Bowman  (Ref.  52)  to  include 
reactions  9,  10,  11,  29,  26,  20,  21,  and  4.  Similar  mechanisms  have  been  proposed  by 
D'Souza  (Ref.  73)  and  Seery  (Ref.  74)  with  D'Souza  also  including  reactions  30,  37, 
and 
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45  HO2  +  HO2  =  H2O2  +  O2 

46  H202  +  H  =  H20  +  OH 


Species  HO2  and  H2O2  are  considered  important  in  low  temperature,  high  pressure,  lean 
flames  (Refs.  75  through  78).  A  common  point  of  agreement  in  Ref.  52  and  Ref.  66 
seems  to  be  that  a  global  hydrocarbon  reaction  (Ref.  79)  is  adequate  to  represent  pyrolysis 
of  hydrocarbons  for  heat  release  rate  calculations: 

47  C„  Hm  +  n/2  02  ->  nCO  +  m/2  H2 

Prediction  of  NOx  concentrations  has  been  shown  to  be  excellent  for  premixed  or 
flat  flames  with  hydrogen,  carbon  monoxide,  methane,  and  propane  fuels,  using  reactions 
2  and  3  (Refs.  66,  72,  73,  74,  80,  and  81).  An  important  point  in  NO  production  is 
the  apparent  high  conversion  of  fuel  N  (as  pyridines,  quinolines,  pyrroles,  indoles  or 
carbazoles)  as  shown  in  Refs.  69  and  82.  Stirred  reactor  experiments  (with  separate  fuel 
injection)  have  shown  disagreements  of  up  to  an  order  of  magnitude  (Figs.  13  and  14) 
between  measured  and  calculated  NO  levels  with  methane  (Ref.  52)  and  propane  (Ref. 
68).  The  experiments  with  methane  also  involved  the  effect  of  poor  mixing.  Excellent 
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Percent  Stoichiometric  Air 


Fig.  14  Nitric  Oxide  as  a  Function 
of  Percent  Stoichiometric 
Air/Propane 

agreement  (Figs.  15  and  16)  was  obtained 
for  carbon  monoxide  and  hydrogen  fuels 
(Ref.  68)  and  where  the  fuel  and  air  were 
premixed  (Ref.  41).  If  poor  agreement  is 
obtained  in  a  controlled  experiment  such 
as  with  stirred  reactors,  it  is  not  surprising 
that  less  than  perfect  agreement  is  obtained 
between  engine  emission  data  (Refs.  7,  1 1, 
12,  58,  83,  through  93)  and  predictions 
based  on  the  single  stirred  reactor  model 
of  the  primary  combustor  in  Ref.  57  (data 
from  Refs.  85,  91,  and  92).  The  recent 
results  obtained  with  the  YJ-93  turbojet 
engine  (Refs.  12  and  93),  in  which 
substantially  different  NO  concentrations 
were  measured  by  physically  sampling 
according  to  SAE  A.R.P.  1256  specification 
(Ref.  12)  and  in  situ  absorption 
spectroscopy  (Ref.  93),  are  significant  when 
comparing  predictions  with  engine  exhaust 
plane  data.  To  date  this  difference  has  not 
been  resolved.  It  is  surprising  that  the 


Percent  Stoichiometric  Air 


Fig.  15  Nitric  Oxide  as  a  Function 
of  Percent  Stoichiometric 
Air/Carbon  Monoxide 


50  100  150 


Percent  Stoichiometric  Air 
Fig.  16  Nitric  Oxide  as  a  Function 
of  Percent  Stoichiometric 
Air/Hydrogen 
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models  of  Ref.  54,  62,  and  69  do  as  well  as  they  do  (less  than  an  order  of  magnitude 
disagreement)  since  they  assume  equilibrium  combustion  of  hydrocarbon  and  air  in  eddies, 
well  mixed  on  a  macroscale  but  not  on  a  microscale,  with  assumed  averaging  techniques 
for  NO  production  in  the  eddies.  Equilibrium  combustion  ignores  the  "super-equilibrium" 
concentration  of  O  atoms  predicted  by  finite  rate  models,  which  can  dominate  the  NO 
production  (Ref.  72).  - 

Reactions  not  considered  previously  may  have  a  significance  yet  to  be  demonstrated. 
Hydrocarbons  do  not  pyrolyze  according  to  reaction  47,  but  rather  through  a  complex 
chain  of  hydrocarbon  fragments,  each  of  which  can  reactor  pyrolyze  further  (Refs.  56 
and  94).  A  possible  fate  of  H  atoms  and  hydrocarbon  radicals  (R)  is: 


48  II  +  NO  +  M  -»■  UNO  +  M 

49  HNO  +  H  -*•  H2NO 

50  M  +  HNO  +  HNO  -♦  (HNO)2  +  M  -  H20  +  N20  +  M 

51  R  +  NO  +  M  -*•  RNO  +  M 

52  M  +  R  +  RNO  -*•  R2NO  +  M 

53  R  +  R2NO  +  M  -*  R2NOR  +  M 


NO  has  long  been  employed  as  a  scavenger  of  hydrocarbon  radicals  although  usually  at 
room  temperatures  (Ref.  95).  This  scavenging  ability  may  have  considerable  importance 
in  afterburners.  The  results  of  Ref.  1 2  indicate  that  NO  virtually  disappears  on  the 
centerline  (total  NOx  remained  constant)  of  the  General  Electric  (GE)  YJ93  afterburner, 
where  there  was  also  a  relatively  high  hydrocarbon  concentration. 

No  concern  has  been  given  to  the  important  effect  (Refs.  96  and  97)  of  third  body 
efficiencies  in  the  above  analyses.  Die  majority  of  reactions  are  affected  by  the  presence 
of  other  species  cither  directly  as  in  three-body  reactions  or  indirectly  through  vibrational 
relaxation  of  species.  H20  and  C02,  for  example,  are  efficient  third  bodies  and,  together 
with  traces  of  alkali  metals  (always  present  as  impurities),  rapidly  relax  excited  nitrogen 
molecules  (Ref.  98).  Reaction  14  is  discussed  in  Ref.  99.  Engelmann  (Ref.  68)  has  shown 
that  an  important  source  of  NO  under  fuel  rich  conditions  is  reaction  4,  generating  up 
to  20  percent  of  the  total  calculated  value, 

N  +  OH  =  H  +  NO 

An  alternate  reaction  may  be  considered  (Ref.  96) 

48  H  +  NO  +  M  -»■  HNO  +  M 
which  is  rapidly  followed  by  (Ref.  100) 

54  H  +  HNO  -*  H2  +  NO 
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or,  if  UNO  »  H,  reaction  50.  Since  the  UNO  is  an  excited  state  complex  on  formation, 
vibrational  de-excitation  of  the  molecule  (eliminating  the  reverse  of  reaction  48  - 
dissociation)  could  play  an  important  part  in  the  net  rate  of  reaction  54.  A  similar  alternate 
reaction  is  (Ref.  96): 


44  H  +  CO  +  M  -*■  HCO  +  M 

55  H  +  HCO  H2  +  CO 

Both  these  reactions  and  reaction  30  can  remove  H  atoms  (or  hydrocarbon  fragments) 
from  the  reaction  chain  with  oxygen  molecules  and  atoms  and  thus  could  alter  the  net 
rates  of  the  nitrogen/oxygen  reactions  (more  O  atom  available  in  reaction  2).  Since  some 
gas  turbine  combustors  work  at  high  pressures  (up  to  30  atm),  these  three-body  reactions 
are  favored. 

3.2  SOOT  AND  SMOKE  PRODUCTION 

It  is  necessary  to  examine  the  process  of  fuel  injection  and  evaporation  in  a  turbine 
combustor  to  understand  the  mechanism  for  soot  (smoke)  production  and  consumption. 
The  fuel  is  injected  as  a  liquid  cone  via  high  pressure  atomization  nozzles  producing  sprays 
with  typical  droplet  sizes  ranging  from  10  to  150  n  with  Sauter  mean  diameters  of  about 
70  /x.  The  cone,  which  may  be  solid  or  hollow,  is  generally  swirled,  and  has  an  included 
angle  between  60  and  80  deg.  The  formation  of  a  stable  spray  from  the  liquid  "sheet" 
occurs  within  a  fuel  nozzle  diameter,  and  thus  interaction  with  primary  zone  gases  is 
usually  between  discrete  droplets.  The  liquid  drops  evaporate  at  a  finite  rate  (providing 
the  combustor  pressure  is  less  than  the  critical  pressure  of  the  fuel);  hence,  following 
the  models  of  Refs.  21,  65,  and  101,  part  of  the  fuel  remains  liquid  until  it  reaches 
the  second  reactor  and  recirculation  flow,  i.e.,  it  traverses  the  first  reactor  intact.  All 
models  consider  each  reactor  to  be  essentially  stoichiometric  (and  thus  generate  high  flame 
temperatures)  but  neglect  the  liquid  fuel  fraction  implied  in  models  including  more  than 
one  reactor.  Thus,  if  these  models  are  correct,  ideal  conditions  exist  for  thermal 
decomposition  of  fuel  to  soot  and  hydrocarbon  fragments,  both  in  the  vicinity  of  the 
injector  and  in  the  first  reactor.  Since  only  about  90  percent  of  the  fuel  is  burnt  in 
the  primary  combustor  and  9+  percent  in  the  secondary  (Ref.  102),  some  of  the  fuel 
goes  right  through  to  the  turbine  and  exhaust.  The  effects  of  engine  operating  variables 
such  as  pressure,  temperature,  mass  flow  rates,  etc.,  on  fuel  droplets  and  drop 
lifetimes/residence  times  is  discussed  in  detail  in  Refs.  101  and  103.  The  principal 
conclusion  of  Ref.  103  concerning  smoke  production  is  that  soot  and  smoke  are  primarily 
a  result  of  inadequate  mixing. 

Expressions  for  soot  and  smoke  production  rates  are  not  available  in  terms  of  detailed 
reactions  and  rates.  It  is  possible  to  estimate  these  using  experimental  engine  data,  if  the 
correlations  of  Refs.  8  and  9  are  accepted.  Thus  if  a  relationship  could  be  found  which 
relates  soot  particulates  (measured  as  a  smoke  number)  to  the  hydrocarbons,  CxHy 
(measured  as  equivalent  methane,  CH4)  and  the  relationship  between  CO  and  CxHy  of 
Refs.  8  or  9  is  used,  then  an  approximate  method  of  accounting  for  fuel  fragments  and 
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their  subsequent  chemistry  downstream  of  the  primary  zone  is  available.  Of  course,  no 
information  is  furnished  on  hydrocarbon  fragment  chemistry  in  the  primary  zone.  A 
different  correlation  is  required  for  different  engines  as  discussed  previously. 

Expressions  for  soot  consumption  have  been  determined  (Refs.  104  and  105).  Soot 
consumption  was  correlated  by  the  partial  pressure  of  hydroxyl  radicals  in  Ref.  105. 
Hey  wood  (Ref.  61)  showed  that  the  data  of  Ref.  105  were  taken  over  a  broader  range 
of  oxygen  pressures  but  that  neither  were  in  the  range  of  turbine  operating  conditions. 
This  is  an  area  in  which  research  is  necessary  before  reactions  and  rates  can  be  prescribed 
to  properly  account  (quantitatively)  for  soot  and  smoke  production  and  consumption. 

SECTION  IV 

DISCUSSION 

The  preceding  sections  outline  the  current  state-of-the-art  in  turbojet  combustor 
modelling.  Models  which  comprise  a  series  of  reactors  which  include  combinations  of  stirred 
reactors  linked  with  plug-flow  zones  are  recommended  because  they  alone  arc  capable 
of  predicting  gross  operating  characteristics  correctly.  A  method  of  estimating  the  individual 
reactor  volumes,  volume  flow  rates,  and  particle  residence  times  is  suggested  using  jet 
entrainment  theory.  Unfortunately,  no  recommendations  are  established  for  successful 
emission  predictions,  although  the  list  of  chemical  reactions  is  long.  The  major  problem 
lies  in  the  absence  of  detailed  calculations  which  assess  the  importance  of  (1)  including 
finite  rate  fuel  evaporation  to  a  global  hydrocarbon  combustion'  model,  (2)  including 
reactions  between  hydrocarbon  (fragments)  and  nitric  oxide,  and  (3)  ignoring  both  the 
relative  cool  regions  where  liner  cooling  air  is  admitted  and  the  core  flow  which  is  not 
entrained  by  the  secondary  and  dilution  air  jets.  In  addition  to  these  problems, 
consideration  must  also  be  given  to  the  fact  that  both  temporal  and  spatial  inhomogenieties 
exist  -  as  witnessed  by  the  soot  and  smoke  emission  -from  engines.  Methods  have  been 
developed  to  predict  the  effects  of  incomplete  mixing;  however,  no  multiple  reactor 
calculations  have  been  carried  out  with  any  of  these  methods,  and  inadequacies  have  been 
reported  for  single  reactor  calculations.  An  alternate  approach  is  outlined;  however, 
significant  development  work  will  be  required  before  it  can  be  evaluated.  Finally,  it  seems 
that  no  scheme  is  available  for  detailing  the  process  by  which  hydrocarbons  are  thermally 
degraded  to  reacting  species,  soot  and  smoke.  It  is  evident  that  much  research  needs  to 
be  done  before  accurate  models  of  turbojet  combustor  (and  afterburner)  processes  are 
developed.  The  author  firmly  believes  that  a  major  fault  of  most  models  outlined  previously 
is  the  inability  to  predict  operating  characteristics  such  as  combustor  hysteresis,  particle 
residence  times,  and  the  effects  of  incomplete  mixing.  It  is  suggested  that  the  two  former 
parameters  should  be  attended  to  and  substantial  effort  initiated  to  solve  the  latter  problem. 

SECTION  V 

CONCLUDING  REMARKS 

The  turbojet  combustor  is  a  very  complicated  component,  and  quantitative 
calculations  of  performance  characteristics  and  emission  levels  require  a  complex  analytical 
model.  The  Swithenbank  model  will  generate  the  right  qualitative  results  and  it  remains 
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to  be  seen  if  quantitative  accuracy  is  achieved.  Inclusion  of  incomplete  mixing  is  essential 
if  the  range  of  throttle  settings  from  idle  to  takeoff  is  to  be  included.  Much  work  is 
necessary  on  soot  and  unburned  hydrocarbon  production  rates;  however,  the  method 
suggested  in  this  report  should  provide  a  reasonable  estimate.  If  the  model  (chemical  and 
fluidynamics)  is  shown  to  be  correct  at  any  one  throttle  setting,  then  at  least  some 
confidence  in  the  chemical  model  will  have  accrued.  Without  this  confirmation,  there  will 
be  no  mutually  satisfactory  basis  for  determining  whether  spectroscopic  or  physical  gas 
analysis  gives  the  right  results  when  substantial  differences  are  observed  (as  reported  in 
Refs.  12  and  93).  These  differences  are  important  and  should  be  resolved  since  they  may 
incorrectly  bias  the  results  of  Environmental  Impact  Statements  and  the  Climatic 
Impactment  Assessment  Program. 
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